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Abstract 

In this paper, a cooperative transmission design for a general multi-node half-duplex wireless 
relay network is presented. It is assumed that the nodes operate in half-duplex mode and that channel 
information is available at the nodes. The proposed design involves solving a convex flow optimization 
problem on a graph that models the relay network. A much simpler generalized-link selection protocol 
based on the above design is also presented. Both the proposed flow-optimized protocol and the 
generalized-link selection protocol are shown to achieve the optimal diversity-multiplexing tradeoff 
(DMT) for the relay network. Moreover, simulation results are presented to quantify the gap between 
the performances of the proposed protocols and that of a max-flow-min-cut type bound, in terms of 
outage probability. 

I. Introduction 

A wireless relay network is one in which a set of relay nodes assist a source node transmit 
information to a destination node. Practically the wireless nodes can only support half-duplex 
communication [1], i.e., no nodes can receive and transmit information simultaneously on the 
same frequency band. Different cooperative transmission schemes for systems with half-duplex 
nodes have been proposed in the literature. Fundamentally, these schemes consist of two basic 
steps. First, the source transmits to the destination, and the relay listens and "captures" [2] the 
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transmission from the source at the same time. Next, the relays send processed source information 
to the destination while the source may still transmit to the destination directly. Variants of these 
techniques have been proposed and have been shown to yield good performance under different 
circumstances [1], [3], [4]. 

Assuming channel state information (CSI) at the nodes, an opportunistic decode-and-forward 
(DF) protocol for half-duplex relay channels is proposed in [5]. In [6], the authors present routing 
algorithms to optimize the rate from a source to a destination, based on the DF technique that 
uses regular block Markov encoding and windowed decoding [7], [8], for the Gaussian full- 
duplex multiple-relay channel. The achievable rate of [7] for the Gaussian physically degraded 
full-duplex multi-relay channel has been established as the capacity of this channel in [9]. In 
[10], it is shown that the cut-set bound on the capacity of the Gaussian single source-multiple 
relay-single destination mesh network can be achieved using the compress-and-forward (CF) 
method, as the relay powers go to infinity. 

Some simpler cooperative diversity methods based on network path selection have been 
recently reported [11], [12]. These selection methods include: (i) the max-min selection method 
[11], wherein the relay node with the maximum of the minimum of the source-relay and relay- 
destination channel gains is selected; (ii) the harmonic mean selection method [11], wherein the 
relay node with the highest harmonic mean of the source-relay and relay-destination channel 
gains is selected; and (iii) the selection scheme of [12], in which the relay that can correctly 
decode the information from the source and has the best relay-destination channel is selected. 
These methods achieve a DMT of d{r) = {N — — 2r) for an N node relay network and 
multiplexing gain < r < 0.5. This is close to what the distributed space-time coding protocol 
[13] achieves, when N is large. Unfortunately, these network selection protocols perform poorly 
in high-rate scenarios (r > 0.5). 

We have proposed a cooperative diversity design based on a flow optimization approach for 
a three-node network in [14]. In this design, the source node broadcasts two distinct flows to 
the destination and the relay node respectively during the relay's listen period. Then the relay 
forwards this information using the DF approach while the source may also send another flow 
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of information to the destination during the relay's transmit period. This scheme is shown to 
achieve the optimal diversity order for the three-node relay channel and yield performance very 
close to optimal full-duplex relaying in both low- and high-rate situations. 

Here, we apply this cooperative transmission design to a general relay network, wherein 
wireless links are present between each pair of nodes in the network. As in [14], assuming 
CSI is available at all nodes we use broadcasting (BC), multiple access (MA) and time sharing 
(TS) techniques to formulate a flow theoretic convex optimization problem based on the channel 
conditions. Instead of considering a total power constraint for all the transmitting nodes as in [14], 
we subject each node to a maximum transmit power constraint. This yields a more reasonable 
system model for a general wireless relay network, especially when the number of nodes in the 
relay network is large. The resulting relaying protocol will be referred to as the flow-optimized 
(FO) protocol. To obtain a more practical cooperative design we develop a generalized-link 
selection (GLS) protocol, in which we select the best relay node out of the available ones to 
form an equivalent three-node relay network to transmit the information from the source to the 
destination. The benefit of this, over other network path selection strategies, becomes evident 
when the rate requirement is high. It is shown that the simple GLS protocol is optimal in terms 
of the DMT [15] and yields acceptable performance even when the rate requirement is high. 

Recently, in [16], the authors have shown that compress-and-forward (CF) relaying achieves 
the optimal DMT for the three-node, half-duplex network, and that DF relaying can achieve 
the optimal DMT of the four-node full-duplex network. In this work, we show that the optimal 
DMT can be achieved for a general A^-node (N > 3) half-duplex network using the FO or 
GLS protocols. Here, it should be clarified that we consider that the wireless links between 
each node-pair experience independent Rayleigh fading, and this corresponds to the definition 
of non-clustered networks in [16]. The performances of the FO and GLS protocols are evaluated 
numerically in terms of their outage probabilities for four- and five-node relay networks for 
uniform and non-uniform average power gains. The numerical results motivate the use of the 
GLS protocol for situations where computation complexity is an issue and show a remarkable 
improvement over the max-min selection method of [11]. The proposed designs, based on BC 
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and MA alone, are sub-optimal in general. For a fair appraisal of the proposed protocols, we 
compare the proposed protocols to an upper bound on the maximum rate, derived using the 
max-flow-min-cut theorem [17, Thm. 14.10.1]. 

II. A General Design Using A Flow-theoretic Approach 

We consider an A^-node wireless relay network with a link joining each pair of nodes. Each 
such wireless link is described by a bandpass Gaussian channel with bandwidth W and one-sided 
noise spectral density Nq. We denote the power gain of the link from node i to node j as Zij. The 
link power gains are assumed to be independent and identically distributed (i.i.d.) exponential 
random variables with unit mean. This corresponds to the case of independent Rayleigh fading 
channels with unit average power gains. Moreover, we assume that each node has a maximum 
power limit of P and can only support half-duplex transmission. Note that this model can be 
easily generalized to the case where channels may have non-uniform average power gains (for 
which numerical examples are presented in Section HVT) . and where different nodes may have 
different maximum power constraints. More specifically, the latter case can be converted into 
the uniform maximum power constraints case by absorbing the non-uniformity in the transmit 
powers into the average power gains of the corresponding links. In the sequel, we characterize 
the system in terms of the transmit signal-to-noise ratio (SNR), 5* = J^^F, at the input of the 
links. Time is divided into unit intervals, and BC and MA are applied with a TS strategy that is 
optimized to maximize the spectral efficiency (which we call "rate" hereafter for convenience). 
To avoid interference between concurrent transmissions, a time interval is divided into slots: 

• During the first slot, the source may BC to all the other nodes in the network. 

• During the subsequent slots, a relay may BC to all other nodes (except the source node), 
or it may receive flows from all other nodes (except from the destination) through MA. 

• During the very last slot, the source and the relays may send information flows to the 
destination using MA. 

Note that the forwarding of information by the relays is based on the DF approach. For practi- 
cality consideration, it is assumed that the phases of the simultaneously transmitted signals from 
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different nodes are not synchronized. In general, for the above transmission protocol, there would 
be a maximum of 2(A^ — 2) + 2 = 2N — 2 time slots of lengths ti,t2, - ■ ■ , t2N-2 respectively. 

Next, we describe the optimization problem using a graph-theoretic formulation. Define a 
graph G = {V,E), where V is the set of nodes, E is the set of all links joining the nodes in 
the graph, and associate the vector r to represent the flow rates associated with each link in 
E. Thus, the number of elements in r equals the cardinality of E. For convenience, we write 
G = {V,E,-r). Now denote the source by S, the destination by V, and the relay nodes by 
TZi, . . . ,TZn^2- The slotting of a unit time interval, as described above, yields simpler graphs 
for each time slot, that we call basic graphs. A basic graph is either one in which a particular 
node may BC to several nodes, or in which several nodes transmit via MA to a particular 
node. Thus for a basic graph, we need to include only the links between the nodes that may 
participate during the concerned time slot. For example, assume that the relay TZi broadcasts 
to all nodes other than the source, during the i-th time slot. The basic graph is given by 



Gi = {V,Ei,r^ where V = {<S,7^l,■■■ ,7^^_2,^5}, E, = {7^l7^2, ■ ■ ■ , 7^l7^^_2, 7^lI)}, u = 



In general, the proposed design involves TS between the basic graphs to yield the following 
equivalent graph G corresponding to a unit interval (see [18] for a similar idea): 



where the number of elements in each vector r^ is extended to |Ui-E'i| by inserting zeros 
appropriately. The second equality in ([T]) implies that G can be viewed as a linear combination 
of the basic graphs GjS, with the equivalent set of edges given by the union of the sets Ej, 
and the equivalent flow rate vector given by the linear combination of the individual flow rate 
vectors r^. 

To maximize the data rate from the source to the destination through the relay network, we 
need to consider each cut that partitions V into sets and with S E V"^ and V E V^. 
Clearly, there can be 2^^^ such possible cuts for the A^-node relay network. Let these cuts and 




, where is the flow from node A to node B during the i-th 



time slot. 
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the corresponding cut sets be denoted by Ck, V^, and Vj}, respectively, for A; = 1, 2, ■ ■ ■ , 2^~^. 
Further, for the graph G, for any two nodes A E Vjf and B E Vjf, there exists a cut edge 
AB that crosses the cut. Denote the total flow through cut edge AB in a unit time interval 
by xab = Z]i=r^^AB- Now recall from network flow theory [19] that the maximum flow rate 
from the source to the destination is specified by the minimal cut of the equivalent graph ©. 
Consequently, we arrive at the following convex flow optimization problem that can be solved 
using standard optimization techniques: 

/ \ 



max mm 



\A€Vf,BeVf A(iVi,B&Vi A^y;N-2^B(iV^N-2 J 



(2) 



over all flow allocations x\q and all time slot lengths tj, subject to 

• the non-negativity constraints: x\q, U > for all cut edges AB and i = 1, 2, ■ ■ ■ , 2N — 2, 

• the total-time constraint: ti + . . . + t2N-2 = 1> 

• the power (capacity) constraints: 

- for a BC slot the flow rates should lie in the capacity region of the BC channel with 
the transmitting node having a power constraint of P, 

- for an MA slot the flow rates should lie in the capacity region of the MA channel with 
a maximum power constraint P for each transmitting node, 

• the^ow constraints: considering steady state operation, the total information flow out of a 
relay should equal the flow into the relay in each unit time interval. 

Note that the dependence of the objective function on the channel gains and the time slot lengths 
is implicitly expressed through the capacity constraints. Denote the cut separating S from all the 
other nodes and the cut separating T) from all nodes as Cs and C©, respectively. Then we observe 
that the cost function in Q above can be further simplified to max min {x{Cs), x{Ct>)}, where 

Ar-2 Af-2 

x{Cs) = xsv + Y ^sn, and x{Cv) = xsv + Y ^^^.^ (3) 
j=i i=i 

are the total flows across the above-mentioned cuts Cs and C-d, respectively. To see this, con- 
sider the cut C with = {S.TZi,--- ,TZi}, and = {T^^+i, ■ ■ ■ ,71^^2,1^} for some / E 
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{1, 2, ■ ■ ■ , — 2}. The total flow across this cut is given by 

N-2 I ( N-2 \ 

x{C)=Xsv+ XsTlj + ^T^r'D + Y ^T^^T^j ' 

j=l+l i=l \ j=l+i I 

Now, consider node z for z G {1, 2, ■ ■ ■ , /}. According to the flow constraint for node i, 

N-2 I N-2 I 

xn,v+ Yl ^n,n,+ Y ^11,11^= xsn,+ Y ^n,n,+ Y ^n^n,- (5) 

j=l+l k=l,k^i j=l+l k=l,k=ii 

Summing ^ over all z G {1, 2, ■ ■ ■ , Z} we get 

I I N-2 \ / / N-2 \ 

Y + Y ^ ^ ^^'^^ + ^ ■ 

i=l \ j=l+l I i=\ \ j=l+l J 

Since E!=i a^T^^Te, > 0, combining ©, & and © gives x{C) > x{Cs). Similarly, we 

have x{C) > x{Cx>)- Thus the cost function in ^ reduces to the above-mentioned form. 

III. Generalized-link Selection and Its Optimality 

In this section, we present the GLS protocol and establish the optimality of the FO and 
GLS protocols in terms of the DMT. This is accomplished in three steps. First, we apply the 
FO protocol to the three-node relay network. Next, we propose the GLS protocol based on a 
selection strategy that is sub-optimal to the FO protocol of Section UIl Finally, the optimality of 
the GLS protocol, and thereby, that of the FO protocol, is established. 

A. The Three-node Relay Network 

The three-node relay network consists of a source {S), a relay (7^), and a destination (V). We 
specialize the general design described in the previous section to this three-node relay network. 
A unit time interval is divided into two time slots of lengths ti and t2 with ti + ^2 = 1- During 
the first time slot, S sends (via BC) two flows of rates x^^/ti = xi/ti and x\^/ti = X2/ti 
to V and TZ, respectively, resulting in the basic graph Gi. During the second time slot, TZ and 
S send (via MA) two flows of rates x^^^/'^z = 2^4/^2 and 2;|^/t2 = 3^3/^2 to V, respectively, 
resulting in the basic graph G2- Combining the two basic graphs yields the equivalent graph 



as G = tiGi + ^2^2- Note that the information flow of rate Xi/t2 sent by IZ during the MA 
time slot is from the flow of rate X2/ti it received during the BC time slot. Thus, we have the 
flow constraint The rate for this network is specified by the min-cut which is clearly 

min{(a;i + X2 + 0:3), {xy + X4 + x^)}. Hence, the flow optimization problem is given by: 

max min {(xi + X2 + Xs), (xi + X4 + X3)} (7) 

over flow allocations Xi, X2, X3, X4, and time slot lengths ti, t2, subject to 

• non-negativity constraints: Xi,X2,X3,X4 > 0, ^1,^2 > 0, 

• total-time constraint: ti + t2 = 1, 

• power constraints: Sbc ^ S, xi < tiC{ZsvS), X2 < tiC^ZsnS) for the BC slot, 
X3 < t2C{ZsvS), X4 < t2C{ZnvS), X3 + X4 < t2C{ZsvS + ZnvS) for the MA slot, 

• flow constraint: X2 = X4, 

where C{x) = log(l + x), and Sbc, the minimum SNR required for the source to broadcast 
at rates xi/ti and X2/ti to the destination and the relay, respectively, in the first time slot with 
< ti < 1, is given by (see [14, Lemma 3.1] for proof) 



S 



BC 



_i_(e-iAi - 1) + _l_e-i/ti(e-V*i _ 1) for Zsn > Zsv, 
1) + 7i-e-^/*i(e^'i/*i - 1) for Zsn < Zsv 



For ti = 0, Sbc = 0. Note that for the BC slot, the last two constraints are redundant when 
ti > 0, and complements the first constraint when ti = 0. 

The solution of this flow optimization problem is given in Appendix |Al As mentioned in 
Section H the above optimization problem formulation is different from that in [14] wherein the 
sum of the source and relay powers, required to achieve a certain data rate, is minimized. More 
specifically, when considering individual power constraints for each node, we cannot use part 2 
of [14, Lemma 3.1] to describe the power constraints for the MA slot. This is because doing 
so would restrict the flows X2 and X3 such that the sum of powers expended at S and TZ is 
minimized. On the other hand, in the present problem, the power constraints only dictate that 
the flow-rates should lie in the MA capacity region specified by the maximum power available 



at each transmitting node, for the particular fading state. With this modification in the constraint 
for the MA slot, the solution approach to the above problem needs to be markedly different 
from that in [14] as shown in Appendix lAl The maximum information rate from the source S 
to the destination V for different cases is summarized below: 

a) Zs. > Zs.: The maximum ra,e is X(S) ^ C(Z..S) wid, di.cC Cransmission from S Co V. 



*2/*l 



b) Zsv < Zsn- The maximum rate is X{S) = maxo<t2<t2„,„ log 

1 I ^S-D I ij 

hC{ZsvS + ZnvS) with tamax = C{ZsnS)/ [C{ZsnS) + C{ZnvS + ZsvS) - C{ZsvS)] and 



+ 



Thus for a given power limit (i.e. a given S) at the nodes, relaying is advantageous only when 
Zsv < Zsn- Further, the optimal solution always allocates a non-zero flow to the direct link. 
Also, the relay-destination link gain Znx) does not influence the strategy of transmission (i.e. 
whether to use only the direct link or both the relay and direct links), but only the amount of 
information through the relay link. 

B. Generalized-link Selection 

For the general A^-node relay network, the flow optimization solution can be computationally 
demanding even for moderate values of N . The GLS protocol described below provides a simple 
sub-optimal design to address this complexity issue. In essence, the GLS protocol identifies the 
best relay path out of the possible N — 2 relay paths and considers only the chosen relay along 
with the source and destination to form a three-node relay network, which we call a generalized- 
link from the source to the destination, for information transmission. In other words, the aim is 
to choose the best relay such that the equivalent three-node relay network obtained (containing 
the source, destination and the chosen relay) gives the maximum rate over all possible equivalent 
three-node networks containing the source and destination. More precisely, we need to consider 
the following possibilities: 

1) ZsD > Zsn,, for alH G / = {1, 2, ■ ■ ■ , A^— 2}: From the results of the optimization problem 

([7]), it is clear that the maximum rate would be C{ZsvS) with direct transmission of all 

data from the source to the destination without using any relay. 
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2) There exists a. k E I such that Z^TZk > ^sv'- Let the set of all such node indices be K 
and for all i E I \ K, Zsv > ZsHi- For this case, choose the node IZ'^, as the relay such 
that k' = argmaxjtgi^Xfc(S'), where Xk{S) is the maximum rate for the three-node relay 
network with the source S, the relay TZ^ and destination V. 
In terms of the worst-case computational complexities for the FO and GLS protocols, it can 
be seen that, for an A^-node relay network with > 3, the FO protocol involves a max-min 
optimization over 2(N'^ — 2N + 2) variables (all possible flows and time slot lengths), subject 
to — 1 non-linear and 2(A^^ — + 1) linear constraints, whereas the GLS protocol involves a 
maximum of A^ — 2 maximizations of a non-linear concave function over two variables, subject 
to two linear constraints, followed by finding the maximum of A^ — 2 real numbers with a 
worst-case complexity of 0(A^ — 2). Moreover, for A^ > 3, for the FO protocol, the BC slots 
potentially involve (A^— 1)- and (A^ — 2)-level superposition coding (SPC) or dirty paper coding 
(DPC) implementations for S and the relays respectively, while the MA slots at the relays and 
V may involve a maximum of (A^ — 3) and (A^ — 2) interference cancelation (IC) operations 
respectively. On the other hand, the GLS protocol involves a maximum of 2-level SPC/DPC and 
one IC operation for the BC and MA slots respectively, for any A^ > 3. 



C. Diversity-multiplexing tradeoff 

As in [15], the multiplexing gain r = limg^^oo where S is the SNR and R{S) is the rate 
at an SNR level of S. Following [15], we parameterize the system, in terms of the SNR S and 
the multiplexing gain, < r < 1, with the rate increasing with the SNR as i? = rlog(S'). With 
the parameterization (r, S), the diversity order achieved by the transmission scheme is given by 

,/ X -log Pe{r,S) 

d[r) = hm (8) 

5^00 log b 

where Pe{r, S) is the average probability of error when the SNR is S and multiplexing gain r. 
The following theorem, whose proof is outlined in Appendix E establishes the optimality of the 
the GLS protocol (and hence the FO protocol) in terms of the DMT: 



Theorem 3.1: The GLS and FO protocols achieve the optimal DMT d{r) 
for all < r < 1, for the A^-node half-duplex wireless relay network. 
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= (A^- 1)(1 -r) 



IV. Numerical Examples 

Using the outage probability as the performance metric, we compare the FO and GLS protocols 
against the max-min selection method of [11], as it provides the best performance amongst 
previously proposed path selection methods, and an outage probability lower bound derived 
using the max-flow-min-cut theorem of [17, Thm. 14.10.1]. 

For the four-node relay network, there can be 6 possible time slots in the FO protocol as 
shown in Fig. [T] To derive an upper bound on the achievable rate (and thereby a lower bound on 
the outage probability), we use max-flow-min-cut type bounds for half-duplex communication. 
There are four possible time slots as shown in Fig. [21 with the first BC slot and the last MA 
slot at the destination same as in the FO protocol, but now, the source and a relay may transmit 
simultaneously to the other relay and the destination during each of the intermediate slots over 
interference channels. We use the max-flow-min-cut theorem to upper bound the maximum 
information flow in these two time slots. 

For the five-node relay network, there can be 8 possible time slots in the FO protocol - four 
BC slots for the source and the three relays to transmit information, and four MA slots for 
the three relays and the destination to receive information respectively. Similar to the four-node 
relay network case, for the max-flow-min-cut bound, there are 8 possible time slots with the 
first BC slot and the last MA slot at the destination being the same as for the FO protocol, and 
multi-source-multi-destination transmissions during the six intermediate slots. 

With the above division of time slots, the formalization of the problem is done as in the 
previous sections, and we use the optimization routine of [20] to obtain the maximum achievable 
rates and upper bounds for different values of required rates. In Figs. [3] and |4l we plot the 
outage probabilities of the various schemes with the required rate R at Ibit/s/Hz and 6bits/s/Hz 
respectively, for the four-node relay network. Figs. [5] and [6] present the same for the five-node 
relay network. When compared to the FO protocol, the GLS protocol suffers a loss of around 
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l.OdB, and around 1.5dB (when R is either Ibit/s/Hz or 6bits/s/Hz), at an outage probability of 
10~^, for the four- and five-node relay networks respectively. On the other hand, the performance 
degradation for the max-min selection method of [1 1], as compared to the FO protocol or even the 
GLS protocol, is more than 12dB at an outage probability of 7.0 x 10^^, when R = 6bits/s/Hz for 
the four-node relay network, and an exactly similar situation can be observed for the five-node 
relay network. Moreover, for the four-node relay network, the FO protocol is within 2.14dB 
(when R = Ibit/s/Hz) to within 7.05dB (when R = 6bits/s/Hz) of the lower bound on the 
outage probability when the outage probability is 10^''. For the five-node relay network, the 
corresponding differences are approximately 3dB and 9.6dB respectively. 

The performances of the different protocols for the four-node relay network with non-uniform 
average power gains are presented in Figs. |7] and [8l and Figs. |9] and \W\ for cases A and B 
respectively, with the average power gains as stated in the figures. In case A, the source-relay 
links are, on average, better than the direct link, and one relay node is, on average, a better 
candidate to forward the information to V. On the other hand, in case B, no one relay has 
very good source-relay and relay-destination links, whereas the inter-relay link is, on average, 
very good. This promotes increased inter-relay flows when using the FO protocol, and thereby 
highlights the limitations of the GLS protocol. The differences between the outage performances 
of the FO and GLS protocols, at an outage probability of 10^^, are L2dB or LOdB, and 2.0dB 
or L3dB (when R = Ibits/s/Hz or R = 6bits/s/Hz), for cases A and B respectively. Thus, the 
gap between the FO and GLS protocols decreases as the required data rate increases. When the 
required rate is high, the coding gain offered by a protocol heavily relies on the efficient use of 
the direct link, and since the usage of the direct link is similar for both the FO and the GLS 
protocols, the performance gap narrows as the required data rate increases. On the other hand, at 
the same outage probability, the difference between the outage performance of the FO protocol 
and the lower bound increases from 1.5dB to 7dB, and from 1.9dB to 6.0dB as the required 
rate increases from Ibit/s/Hz to 6bits/s/Hz, for cases A and B respectively. Overall, these results 
demonstrate trends similar to the uniform average power gain case, and confirm the generality 
of the proposed protocols. 
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V. Conclusions 

We proposed a cooperative transmission design for a general multi-node half-duplex wireless 
relay network. It is based on optimizing information flows, using the basic components of BC and 
MA, to maximize the transmission rate from the source to the destination, subject to maximum 
power constraints at individual nodes. We also proposed the simpler GLS protocol, that combines 
relay selection, and flow optimization for a three-node relay network. These protocols were shown 
to achieve the optimal DMT for a general relay network. Simulation results for the four- and 
five-node relay networks for uniform and non-uniform average power gains demonstrate that the 
performance of the much simpler GLS protocol is close to that of the FO protocol. This suggests 
that the GLS protocol can be used in systems with low-complexity requirements. We also note 
that the proposed FO and GLS protocols can be used in wireless networks with topologies more 
complicated that the wireless relay network considered here. For example, application of similar 
ideas to a parallel relay network in which there is no direct connection between the source and 
the destination is considered in [21]. 

Appendix 

A. Solution to optimization problem (|7]) 

We consider two cases with regard to the link gains: (a) Zgv > Zgn^ and (b) Zsv < Zsn- 
For both cases, we solve the optimization problem in two stages: first, we fix ti, ^2 > such 
that ti + ^2 = 1 and find the optimal flows xi, X2, X3 in terms of ti, t2, and then, find the 
optimal values for ti, t2 to maximize the objective function. 

a) Zsv > Zsn^ To obtain an analytical solution to the optimization problem and better 
insight into the nature of the solution to the flow optimization problem, we modify the repre- 
sentation of the BC slot power constraint from that in ([7]) to the one that is more conventionally 
used to describe the capacity region of the Gaussian BC channel, as presented in Using the 
flow constraint in (|7]), we first solve ^ for fixed ti, ^2- 



max(a;i + X2 + 2:3) over xi, X2, X3, a, subject to 



(9) 
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Xi, X2, Xs > 0, < a < 1, 

ZsnaS 



xi < tiC{ZsvaS), X2 < tiC . , , „ ^ , , 
X2 < t2C{ZnvS), X3 < mill {t2C{ZsvS), t2C{ZsvS + ZnvS) - X2} ■ 

Here, a G [0, 1] is tlie fraction of total power spent at the source to transmit xi directly to the 
destination during the BC slot, and a = 1 — a. Although, this modification of the BC slot power 
constraint apparently makes the optimization problem non-convex owing to the non-convexity 
in a, as we shall see in the sequel, this issue can be handled easily by utilizing the monotonicity 
of the logarithm function. 

Denote the optimal solution by {x\, X2, x*^, a*) and the corresponding maximum rate by 
X{ti,t2). It is clear that x\ = tiC{Zsva*S). Suppose that xl < hC ( i+g^^fg )- Since 
tiC ( i^f^^fg ) is a decreasing function of a, we can increase a from a* to such that 
Xi = tiClZsva'^S) > xl and x^ = tiC ( i^ff^^aog )- Thus the objective function becomes 
{x\ + X2 + a;3) > {x\ + X2 + xl) at This contradicts the optimality of {x\, x*2, x*^, a*) . 
As a consequence, we have x*2 = hC [j^f^^) < t2C{ZTivS). This implies that a* > 

,0 /•, and a;° < 1. In essence, this means that the optimal 



1 



max 



ZsnS 



1 + ZsnS _ ^ 



[I + ZnvSy-l'^ 



Xl and X2 should lie on the boundary of the degraded BC capacity region. With this, it is 
obvious then that S3 = min {^2^(^51,5), t2C{ZsvS + Z^vS) - hC ( i+g^fg )}- Therefore 
the optimization problem of ^ can be re- written as: 

max (xi + X2 + xs) (10) 

subject to max{0, a°} < a < 1, Xi = tiC{ZsvOiS), X2 = hC ( ^^—~~r, 

V 1 + ZsTiOiS / 

X3 = min !^t2C{ZsvS), t2C{ZsvS + ZnvS) - hC (^f^^) 
We observe that X3 = t2C{ZsvS) above if and only if 
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Comparing this to the expression for a[ 



grves a, 



'° < < 1. 



Next we consider two possible sub-cases: 




can be expressed as X{ti,t2) = max{Xi(ti, ^2), -^2(^1,^2)} where 



— hmax- The maximum rate 



^1(^1,^2) 



max tiC{ZsDaS) + t2C{ZsvS + Z-ji-pS) 



= tiC{Zsv<S) + t2C{ZsvS + ZnvS) 



X2{tl^t2) 





where the first inequality in (fT2l) holds since (fTT)) is not satisfied, and the second inequality 
in (fT2l) holds since ol\ < I and that the first term in the previous step is monotonically 
increasing in a\ when Zsv > Zsn- This way, the last observation helps avoid the non- 
convexity issue mentioned before. Similar arguments hold for X2(ti,t2) in (fT3l) . 



Hence, (fT2l)-(fT4l) imply that the solution to (|9l), for any ti, t2 pair, occurs at a = 1, and the solu- 
tion to the original problem of (|7]), when Zsv > Zsn, is given by max{o<ti,t2 ; 1^+12=1} -^(^1, ^2) = 
C{ZsvS) with any ti, ^2 pair such that ti, ^2 > and ^1+^2 = 1- 

Z7j Zsv < ZsTi: In this case, the source-relay link is better than the source-destination link. 
Again, we first fix time slot lengths ti and t2 and solve for the optimal values of xi, X2, X3, and 
a, and then maximize the objective function of © over all feasible time slot lengths. Following 



ii. t2C{-^^)>hCiZsnSy. 



In this case, we have < and t2 > t2max- Thus the maximum rate is given by 
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similar arguments as in Case a), the optimization problem of ([7]) can be re-written as: 



max (xi + X2 + X3) 
subject to < a < min{l, 0;°}, Xi = tiC 



(15) 



ZsvCiS 



X2 = tiC{ZsnaS), 



,1 + ZsvOiS , 

X3 = min{t2C{ZsT:)S), t2C{ZsvS + Z^vS) - tiC{ZsnOiS)} 



^ ^ I+Zs-dS) 



A -1 
= "6- 



where a° = ^^^^ (1 + Z-jidSY^^^'^ — 1 is an upper bound on a such that tiC^ZsnOiS) < 
t2C{Z'jivS). As before, X3 = t2C{ZsvS) if and only if a < 
Also, < < Again, we consider two possible sub-cases: 



In this case we have al < 1 and to < , ^ 



C{Zs-R.S) 



X(ti,t2) can be expressed as X{ti,t2) = raax{Xi{ti,t2), X2{ti,t2)} where 



hmax- The maximum rate 



Xl{h,t2) 

-^2(1^15 ^2) 



max tiC 



0<a<al 



ZsvOtS 

1 + Zsx)aS 

ZsvCiS 



max tiC , _ 

aJ<«<min{l,«5J} V 1 + ^5©Q;o 



) + tiC{ZsnaS) + t2C(Z5755), (16) 
t2C{ZsvS + ZnvS), (17) 



and both the maxima in (fT6l ) and (flTI) are attained at a = a^. Substituting the expression 
for al in (fT6l) and (fTTI) . we obtain Xi(ti,t2) = -^2(^1,^2) and 



/ 



X(ti,t2) =tl log 



1 + ZsvS 



\ 



1 _|_ Z^T, 

Zsn 



1 I Z-!ZT}S 



i2/ti 



^2C(Z5I,5 + Z7^^,5). (18) 



Jn this case we have al > I and ^2 > hm^x- Thus the maximum rate is given by 

X{h,t2) = me^x hc(-^^^^)+hC{ZsnaS)+t2C{ZsvS) 

0<a<l \1 + ZsTidb / 

= hC{ZsnS)+t2C{ZsvS) 



(19) 



where the maximum occurs at a = 1, as the function to be maximized is monotonically 
increasing in a when Z^-jz > Zsv- 
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Finally, we optimize the above solution to (fTSl ) over all possible time slot lengths to obtain the 
solution to the original problem in (|7]) when Zsx> < Zsn- Corresponding to Case i. above, when 
^2 < t2max, we note that 



max XitxM) > X(l - ^2,^2)1*2=42, 

{0<ti,t2:ti+i2=l} 



max 



_ C{ZsnS)C{ZsvS + Z-jiuS) 

C{ZsTiS) + C{Z'jix>S + ZsdS) — C{ZsvS) 

On the other hand, corresponding to Case ii., when t2 > t2max, from (fT9l ), we have 



(20) 



max X{ti,t2) < {I -t2m3x)C{ZsTlS) +t2m^xC{ZsvS) 

{0<ti,t2: ti+t2 = l} 

C {Z snS)C {Z svS + Z-juyS) 



(21) 



C{ZsnS) + C{Ztix)S + ZsdS) — C{ZsvS) 
where the inequality in (|2T)) is obtained from (fT9l) by using Zsn > Zsv and t2 > t2max- Hence, 



from (1201) and (|2T]) . we conclude that when Z^ji > Zsv, the maximum achievable rate is given 

l+Zs-yS 



by X{S) = maxo<t2<t2^,, h log 



1 I £SE. 



t2C{ZsvS + Z'jix>S)- 



B. Proof of Theorem I3.il 

In this appendix, we sketch the proof of the theorem. Part 3) of Theorem 4.2 of [14] can be 
generalized for the A^-node relay network to prove that, as the block length goes to infinity (during 
any particular time interval), the average error probability for the FO protocol is upper bounded 
by its corresponding outage probability. Here, the outage probability denotes the probability 
that the data rate R cannot be supported by the system when the SNR is S, i.e., Poa\.{r, S) = 
Fr[X(S) < rlogS] where X(S) denotes the maximum rate possible for the given channel gain 
realizations when the SNR is S. Thus, from the definition of diversity order ([8]), we have 

,/ X ^ v -logPout(r, S) 

d[r) > hm . (22) 

S^co log S 

Moreover, the above result from [14] can be directly used to prove the same for the GLS protocol. 
Using this fact, we derive a lower bound to the DMT that can be achieved by the GLS protocol. 
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The sets / and K, used in the sequel, are the sets of indices as described in Section IIII-BI The 
outage probability for the GLS protocol is given by 

P!Ur, S) = Pr (m^xXkiS) < r log (23) 

where Xk(S) is the maximum rate achievable by the three-node relay network formed by the 
source S, the relay TZ^ and the destination V. We have the following possibilities: 

• Case A: = 0, i.e. the cardinality of the set K is zero. This corresponds to the case 
when > ^57^fe for all A; G /. 

• Case B: li^l = i with i E I. 

Note that for Case B there are (^7^) possibilities for the set K with cardinality i. Since the 
link gains are assumed to be i.i.d., and the outage probability depends on the distribution of the 
maximum of Xk(S) over all /c G / (or effectively, over all k E K when \K\ > 0), only the 
cardinality of K is significant. Let the [~i possible constructions of the set K be represented 



by a "generic" set Kq with cardinality i. Without loss of generality, we describe Kq as the set 
K corresponding to the case when the indices of the relay nodes are ordered according to their 
source-relay link gains, i.e. Zsn-^ > Zsn^ > ■ ■ ■ > Z 311^-2- Thus, Case B now implies a solitary 
choice for set K, viz. Kq = {1,2, - ■ ■ Therefore, from (|23] ). we have 

P!Ur,S) = FT{C{ZsvS)<rlogS\\Ko\ = 0)FT{\Ko\=0) 

N-2 , 

-|- Pr ( max Xfc (5*) < r log 5* 



K^\=i\Vi{\K^\=i). (24) 



We observe that using the right-most expression of (|20l) instead of Xfc(S'), for each k E Kq, 
in (|24l) gives an upper bound on Poltir,S). This is utilized in obtaining a lower bound on 
the diversity order of the GLS protocol. Let {Sn}'^=i be an increasing unbounded sequence of 
SNRs with 5*1 > 1. Define the sequence of random variables {M^jJ^^, {B^^}^^^ and {A^j^^ 
"i'h = "l^SS:^:^' . Bt = and M = ^ = t^, respectively. 

Note that for all k E Kq, Q a.s. This implies that {B^ - A^) ^0 a.s. Define 

A'^ = mo^KkfzKoA^ and B'^ = max^g/^g 5^ (= 5^). Then using the above, it can be seen 
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that (B'^ — A'^) a.s. Further, lim Pr (B'^ < r\ \Kq\ = i) exists, and therefore the above 
implies that lim Ft (A' < r\ \Ko\ = i) = lim Fi (B'„ < r\ \Kq\ = i). Using this in (l22l) and 
((24)) . the diversity order for the GLS protocol, d^^(r), satisfies 

- log P!Ur,S) 



d^\r) > lim 

5— »oo 



- ^™ ] — \r 

n-»oo log Sr, 



logS 
log 



Pr ( Zsv < 



or 



Kol =0 PrdiTol =0) 



N-2 , — 1 

i=i ^ 
lim - — ^log 

log Sn 

N-2 



\Ko\=t)Fi{\K, 



Sri 



Ft {max{Zsx>,ZsTi^, 



S^ 1 



Sn 



(25) 

li^ol =0)Pr(|i^^o| =0) 



i=l 



or 



J2 v[iax{Zsv, ZsTii, ■ ■ ■ , ZsnN_J < — 



Sn 



\Ko\ = t)Fii\K, 



(26) 



lim 

S^oo 



lim 



- log (Pr (maxjZ^i,, Zsm, ■ ■ ■ , Zshm^i} < 

log 5 

-log([l-exp(-(5^'-l)/5)r-i 



log 5 



(iV- 1)(1 -r) 



(27) 



where (|26|) is obtained from (l25l) by noting that maxjZ^x), ^57^^, • • • , ZsnN-2} = when 
|i^o| = 0, and max{Z5x,, ^57^^, ■ ■ ■ jZ^t^^^J = ma^XkciKo Zsn^ = Zsui when |A"o| > 0, the 
first equality in (l27l) is due to the link gains being i.i.d., and the second equality in (|27l) is 
obtained by using L' Hospital's rule. 

Next given an A^-node relay network, consider the multiple access cut that separates the 
destination from all the other nodes. Clearly, the total flow across this cut gives an upper bound 
on the maximum rate achievable in the A^-node relay network. Consequently, a lower bound on 
the outage probability Po^^^{r, S) can be obtained using the maximum sum-rate across this cut: 



PUr,S) > Ft C{{Zsv + Zn,v + --- + Zn^^^j,)S)<r\ogS 



Pr 



{Zsx> + ^TeiU + 1- Zti^_2X)) < X 



l{N - 1, {S^ - 1)/S) 
T{N-1) 



(28) 



where 7(0, x) = Jq t"' ^dt is the lower incomplete gamma function and T{a) = /q°° ^dt 
is the complete gamma function. The result in part 1) of Theorem 4.2 of [14] can be extended to 
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show that the diversity order of any transmission scheme over the wireless relay network must 
satisfy 



Finally, from (|T7l) and (|29l) . we see that the GLS protocol, and hence the FO protocol, achieve 
the optimal DMT of (N - 1)(1 - r) for all < r < 1. 
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Fig. 2. Transmission strategy to obtain a lower bound on the outage probability for the four-node relay network. Here 
ti + ■ ■ ■ + 14 = 1, and the optimization is over xi, ■ ■ ■ , Xi4, and ti, ■ ■ ■ , ti, with the application of the max-flow-min-cut 
theorem for the intermediate slots. 
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SNR (dB) 



Fig. 3. Four-node relay network with uniform average power gains: Outage probabilities for required rate R — Ibit/s/Hz. 




SNR (dB) 



Fig. 4. Four-node relay network with uniform average power gains: Outage probabilities for required rate R — 6bits/s/Hz. 
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Fig. 5. Five-node relay network witli uniform average power gains: Outage probabilities for required rate R = Ibits/s/Hz. 
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SNR (dB) 



Fig. 6. Five-node relay network with uniform average power gains: Outage probabilities for required rate R = 6bits/s/Hz. 
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Non-uniform average power gains: Case A 




2 4 6 8 10 12 14 
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Fig. 7. Four-node relay network with non-uniform average power gains. Case A: £[^57^^] = 2.0, £[^57^3] = 2.0, E[Z5-p] = 
1.0, ¥,[Z'jZiiZ2] = 1-0, IE[^7Ji'd] = 1.5, ¥,[Ztz2'd] ~ 1-0. In this case, both the source-relay links are, on average, better than 
the direct link, and one relay-destination link is, on average, better than the other, resulting in relay TZi being more favorable 
than relay TZ^- 
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Non-uniform average power gains: Case A 
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Fig. 8. Four-node relay network with non-uniform average power gains. Case A: EfZ^Tj^] — 2.0, £[^57^3] = 2.0, E[Zs-d] = 
1.0, ¥,[Z'jZiiZ2] = 1-0, IE[^7Ji'd] = 1.5, ¥,[Ztz2-d] ~ 1-0. In this case, both the source-relay links are, on average, better than 
the direct link, and one relay-destination link is, on average, better than the other, resulting in relay TZi being more favorable 
than relay TZ2- 
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Non-uniform average power gains: Case B 
Spectral efficiency, R=1 bits/s/IHz 




SNR (dB) 



Fig. 9. Four-node relay network with non-uniform average power gains. Case B: E[^57jJ = 1.5, E[Z5-;j,2] = 0.75, EfZs-p] = 
1.0, E[Z-R.j7j.2] = 3.5, E[Z-r.j-d] = 0.2, ¥,[Z'iz2'd] = 3.0. In this case, one source-relay link is, on average, better than the 
direct link, which, in turn, is better than the other source-relay link. The reverse is true for the relay-destination links, and the 
inter-relay channel is, on average, very good. This situation promotes inter-relay interactions for the FO protocol, and thereby 
highlights the sub-optimality of the GLS protocol. 
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Non-uniform average power gains: Case B 




SNR (dB) 



Fig. 10. Four-node relay network with non-uniform average power gains. Case B; E[Z57jJ = 1.5, £[^57^3] = 0.75, E[Z5'd] — 
1.0, ¥,[Z'iZiiZ2] = 3.5, EIZ-jzi-d] = 0.2, ¥,[Zh2T>] = 3.0. In this case, one source-relay link is, on average, better than the 
direct link, which, in turn, is better than the other source-relay link. The reverse is true for the relay-destination links, and the 
inter-relay channel is, on average, very good. This situation promotes inter-relay interactions for the FO protocol, and thereby 
highlights the sub-optimality of the GLS protocol. 



